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TECHNICAL NOTE 2110 

INTERFEROMETER CORRECTIONS AND MEASUREMENTS OF LAMINAR 
BOUNDARY LAYERS IN SUPERSONIC STREAM 
By Robert E. Blue 


SUMMARY 

The laminar boundary layer on a flat plate in a supersonic 
stream was investigated by means of a Zehnder-Mach interferometer 
and a total-pressure probe. The density gradient in the boundary 
layer on the plate necessitated consideration of the effect of light 
refraction on the interference data. 

Density and velocity distributions, which were obtained with 
the two instruments for a range of local Reynolds number from 
0.318 x 10® to 1.08 X 10® at a nominal Mach number of 2.02, compared 
favorably with laminar-boundary- layer theory. Skin-friction coeffi- 
cients calculated from the boundary- layer profiles were larger than 
predicted by theory. 


INTRODUCTION 

Increasing application of light-interference methods to aero- 
dynamic and thermodynamic studies has indicated that the procedure 
of obtaining and evaluating interference photographs differs con- 
siderably, depending on the study to which the method is applied. 
Factors such as light refraction, definition of a solid surface, and 
nonuniformities in the density field may be of little or no concern 
in one application, for example, when measuring pressure distributions 
on an airfoil in subsonic flow (reference 1); whereas they may cause 
an appreciable error for another application, such as the measurement 
of temperature distribution around a heated body (reference 2) . Even 
when interferometry is applied to the specific problem of obtaining 
density distributions in a boundary layer of air near a solid surface, 
these factors can cause errors of varying degrees depending on the 
type of density profile (laminar or turbulent), the boundary- layer 
thickness, and the total density change across the boundary layer. 
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Boundary- layer investigations on an airfoil and a flat plate 
in subsonic flow were made at the Hermann Goring Institute by Zobel 
in 1940 (reference 3). The boundary layers were relatively thick 
(of the order of 0.1 in.) and of a turbulent nature, with a total 
density change across the layer of approximately 2.8 percent of the 
density at the edge of the boundary layer. Errors in the method 
apparently were assumed to be negligible, inasmuch as no discussion 
was devoted to them. 

A boundary- layer investigation, conducted at the NACA Lewis 
laboratory and presented herein, is similar to the earlier work done 
in Germany, except that errors caused by light refraction and an 
inability to distinguish the exact location of the surface are 
amplified because the boundary layer studied is much thinner 
(0.030 in.), the density change larger (42 percent), and the pro- 
file of a laminar type. A laminar boundary layer concentrates most 
of the total density change across the layer in a relatively narrow 
region, whereas for a turbulent profile the density change is more 
evenly distributed. The high density gradients that exist in a 
laminar layer amplify errors caused by light refraction. 

Although some information was obtained regarding the behavior 
of a laminar boundary layer on a flat plate in a supersonic stream 
both by interferometer and pressure-probe measurements, the experi- 
ments and the analysis were made with the intention of systematically 
organizing the steps required to evaluate and to correct interference 
data and to estimate the error in the final result. With this infor- 
mation, a more detailed study of heat-transfer and boundary- layer 
phenomena can be undertaken using the interferometer as the primary 
measuring instrument. 


APPARATUS AND PROCEDURE 
Interferometer 

A Zehnder-Mach type interferometer was used for this investi- 
gation. The construction and the operation is conventional and will 
not be discussed in detail as there are many excellent descriptions 
of the problems of design, construction, and operation of this type 
of instrument . (For example, see references 4 and 5.) 

A sketch of the interferometer and wind-tunnel installation is 
shown in figure 1. The collimator consists of an f/l.3, l/2-inch 
focal-length lens focusing the light from a high-pressure mercury- 
vapor light on a 0.008- by 0.125-inch slit located at the principal 
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focus of an f/2.5, 12.5-inch focal-length lens. The camera lens is 
an f/7.5 lens with an 18.5-inch focal length. The light is filtered 
for the mercury green line 5461 A. 


Supersonic Wind Tunnel and Model 

The measurements were made on the boundary layer developed on 
a flat steel plate 4 inches long and 3.6 inches wide. The plate had 
a 12° wedge angle on the bottom. The leading edge was honed to give 
the sharpest edge possible. The angle of attack could be varied 
between ±3°. The plate completely spanned the 3.6-inch square test 
section of a two-dimensional supersonic tunnel designed by the method 
of characteristics for a Mach number of 2.08. Flow observation was 
through l/2- inch, optically flat windows set in aluminum side plates. 
The wind tunnel and flat plate with one side plate removed is shown 
in figure 2. Provisions were made for inserting a total-pressure 
tube, which is shown above the plate in figure 2, to probe the bound- 
ary layer. Distance from the surface to the bottom of the probe tip 
was measured with a micrometer and is considered accurate to less 
than 0.0005 inch. Contact of the probe with the surface was indi- 
cated by an oscilloscope, which also provided a method of detecting 
probe oscillation. The oscilloscope was wired in series with a 
battery; one end of the circuit was connected to the plate and the 
other end to the probe. A flip in the oscilloscope trace occurred 
when the plate and the probe, which were electrically insulated from 
each other, came into contact. The pressure measured by the probe 
was assumed to be the pressure at the geometric center of the probe 
opening. The over-all height of the probe tip was 0.0026 inch. The 
shape of a typical probe opening is shown in figure 3. Twelve 
static-pressure taps and three thermocoujfLes in the surfa.ce pro- 
vided a method of determining the nature of the pressure gradient 
and the temperature along the plate. 


Air Supply 

The air supplied to the tunnel had a dew point of 410° E or 
lower. Inlet pressure and temperature could be varied over the 
ranges of 11 to 40 inches of mercury absolute and 540° to 590° E, 
respectively, which gave a stream Eeynolds number range of 1.04 x 10 5 
to 4.27 x 10 5 per inch in the test section. 
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Preliminary Adjustments 

Before any boundary- layer data were taken, a preliminary experi- 
ment was made in which the angle of attack of the plate was varied 
over the range of ±3° in order that a zero pressure gradient would 
exist over the length of the plate. The change in pressure distri- 
bution was negligible over this range of angle of attack. A second 
criterion for adjustment therefore was that the shock wave from the 
leading edge of the plate have a minimum strength. The interferom- 
eter was used to observe the wave while the angle of attack was 
continuously changed over the entire range of possible angles. The 
plate was then set to the angle where the disturbance appeared to 
be the weakest. This adjustment was made when the tunnel was oper- 
ating with inlet air at a pressure of 40-pound gage rather than 
5 -pound- gage used when measurements were made. The purpose of the 
higher pressure was to amplify the density difference and the changes 
in the density difference across the shock wave. Dry 40-pound air 
was generally unavailable for quantitative measurements. An illus- 
trative interf erogram of the leading-edge shock wave is shown in 
figure 4. The strength of the disturbance on the upper side of the 
plate (the side on which the boundary layer was investigated) can 
be compared with the compression region resulting from the turning 
on the underside of the plate. Pressure distributions (fig. 5) were 
recorded at various inlet pressures for the final orientation of the 
plate . 

During the experimental portion of this investigation, it was 
noted that extreme care must be taken to prevent regions of local 
turbulence from developing at the surface of the plate. Such regions 
can be caused by leaky static-pressure orifices or by specks of dirt 
on the surface that are so small as to be almost invisible. 


ANALYSIS 

Total-Pressure Data 

The conventional assumption that the static pressure is constant 
through the boundary layer along a line normal to the surface was 
used in calculating the Mach number distribution from the ratio of 
measured total pressure to wall static pressure by means of the 
Rayleigh equation for Mach numbers greater than 1 and by the 
isentropic-flow equation for Mach numbers less than 1. Mach number, 
ambient temperature, and total temperature are related at any point 
in the flow by the one -dimensional energy equation, from which the 
following expression is derived: 
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_t_ ^ T / 1 + 0.2 Ml 2 \ 

*1 T 1 \1 + 0.2 M 2 J 


(All symbols used in this report are defined in appendix A.) The 
ratio of specific heats was taken as 1.4. In order to calculate 
ambient- temperature distribution from the measured Mach number using 
equation (1), the distribution of the total temperature must be 
known. Laminar-boundary- layer theory has shown that only for the 
case of a Prandtl number Pr equal to 1 can the assumption be made 
of constant total temperature through the boundary layer. It is 
known, however, that Pr ~ 0.72 for air. Inasmuch as no means were 
available to measure the total-temperature distribution, the theory 
developed in reference 6 was used to calculate the total temperature 

/ u.. 

as a function of the dimensionless distance variable y a / — — for 

V U 1 X 

a Mach number M]_ measured at the edge of the boundary layer by the 
total-pressure probe. A typical total-temperature distribution is 
plotted in figure 6. Values of total-temperature ratio taken from 
the theoretical curve and measured values of Mach number were sub- 
stituted in equation (1) to obtain t/t^. The density ratio p/p^ 

is given by 



because p is assumed constant. The velocity ratio u/u ^ was 
calculated from 


u 

U 1 


t 

*1 



(3) 


Interference Data 

The method used to analyze an interferogram is determined by 
the basic adjustment of the instrument. Adjustment to give inter- 
ference fringes parallel to the line along which the density dis- 
tribution is to be measured, in this case perpendicular to the plate 
surface, produces a photograph that requires the least work to ana- 
lyze. A vertical-fringe adjustment was used in this case, inasmuch 
as the flat plate was horizontal. The indicated boundary- layer pro- 
file was calculated by a method described in reference 2. 
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The principle of the method of analysis is illustrated in fig- 
ure 7. If the density were constant throughout the field, the center 
line of fringe A in the figure would pass through point P and would 
be everywhere parallel to the dotted line that is perpendicular to 
the surface. The shift of one fringe at P indicates that there is 
one less wavelength (less because the density decreases) in the path 
of the light appearing to come from P than there would be if the 
density at P were the same as at Q. The difference in density Ap 
between the points P and Q in terms of the fringe shift N (refer- 
ence 2) is 

^0 

40 = E N (4) 

For boundary- layer studies, the density at the point P is associated 
with the distance of the point from the surface. This distance was 
measured directly from the interferogram negative by means of a com- 
parator and was converted to actual distance from the surface by 
dividing by the magnification of the photograph. 

In order to calculate the ratio of boundary- layer density to 
density at the edge of the layer, the absolute density at the edge 
of the layer must be known. This density can be calculated either 
from pressure and temperature measurements or from a measure of the 
number of monochromatic fringes that traverse the position where the 
density profile is to be evaluated when the tunnel changes from a 
no-flow to a flow condition. The method using the monochromatic 
fringes is inaccurate when the boundary layer on the tunnel windows 
occupies an appreciable portion of the tunnel span. The stream 
density was therefore calculated from the wall static pressure and 
the calculated free-stream ambient temperature. The density pro- 
file obtained in the manner just described is hereinafter referred 
to as the "indicated" density profile, as distinguished from the 
actual density profile. 


Errors in Interference Data 

The indicated density profile may or may not be an accurate 
representation of the profile existing in the boundary layer on the 
plate. In order to determine the degree of approximation to the 
actual density distribution, an analysis was made of the factors 
affecting the formation of boundary- layer interference fringes when 
the density profile has a zero density gradient at the surface and 
at the edge of the boundary layer and rises monotonically to a 
maximum value between the zero values. This analysis is perhaps 
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more valid for a laminar profile than for a turbulent profile. The 
analysis indicates the necessary steps required to correct inter- 
ferometer data and also gives an estimate of the error in the final 
results. The errors have been divided into three general catagories: 
those caused by (a) light refraction, (b) interferogram evaluation, 
and (c) tunnel-window boundary layer. 

Eefraction errors. - The principal cause of errors in the 
interference data is the refraction of the light in the nonuniform 
density field of the boundary layer. When the density at the surface 
is less than in the stream, the light is bent away from the surface, 
as schematically shown in figure 8. The path of a light ray through 
the boundary layer can be calculated in terms of the coordinates z 
and y from the law of refraction for non- homogeneous media, which 
was obtained from Fermat* s principle. If density variations are 
assumed to exist only in the plane of the light path, this law can 
be expressed by 


y" 1 dn 1 dn , x 

777? =5 3? - 5 Si y 

where the refractive index of air n is related to air density by 
the experimentally determined law of Biot and Mas cart 



0.1166 (cu ft/slug) 


( 6 ) 


when 


X Q = 5461 A 

If light were not refracted but passed through a two-dimensional 
boundary layer parallel to the surface, and if all other errors were 
negligible, the interferometer would indicate a density profile 
identical to the actual profile. Refraction, however, causes the 
following errors: 

1. Averaging errors. Instead of traversing the boundary layer 
at a constant density, the light passes through layers of increasing 
density. The average density along the light path is therefore 
higher than the density at the entrance height. The difference 
between the local- stream and average densities is therefore less 
than the difference between the local-stream and entrance densities. 
Because the density difference is directly proportional to fringe 
shift (equation (4)), the fringe shift is less than would exist with 
no refraction. 
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2. Window errors. The light, after passing through the bound- 
ary layer, makes an angle 0 e with the plate surface at the inside 
of the tunnel window. On striking the glass at an angle, the light 
is refracted according to Snell's law (sin 6 e = ng/nair sin ®r)> 
and passes through the glass in a straight line at an angle G r 
with the surface. At the outside surface of the window, the light 
is again refracted hack to the original angle G e . Without refrac- 
tion in the boundary layer, the light would travel a distance T 
(the glass thickness) through the glass, or an optical-path length 
of ngT. With refraction, the optical path becomes ngT/cos 0 r , 
which is an increase in optical path over the undeviated path by 
the amount ngT(l - cos 0 r )/cos 0 r . This quantity, divided by 

the wavelength of the light in a vacuum, is the decrease in fringe 
shift between a ray through the boundary layer and a ray through 
the stream. Ths density indicated by the interferometer will there- 
fore be higher than it would be if the effect were not present. 

(The boundary layer and the glass window have an opposite effect 
on optical-path length.) 

3. Displacement errors. As a general rule, the light will not 

appear to come from the point in the object plane at the distance 
from the surface at which the light entered the tunnel. Instead it 
will appear to originate at the point y a (fig. 8) where the back- 
ward extension of the light ray from the outside of the tunnel win- 
dow intersects the object plane of the camera lens. The location 
of this point is influenced by: (a) the emergent height of the 

light y e , (b) the emergent angle 0 e , (c) the thickness of the 

tunnel window T, and (d) the location of the object plane. 

The indicated density is therefore an average density along a 
curved light path and is increased in value by the additional path 
through the tunnel window, and appears to originate at a somewhat 
arbitrary height above the surface, determined by the refraction 
of the light in the boundary layer and the geometry of the optical 
system. 

An estimate of the approximate magnitude and relative importance 
of each of these errors was made by using equation (5) to determine 
the light path through a boundary layer having a density distribution 
corresponding to a cosine curve. The assumed cosine profile resem- 
bles a laminar density profile in that it has the previously men- 
tioned characteristic of zero density gradient at the plate surface 
and at the edge of the boundary layer. The density at the edge of 
the layer, the ratio of stream density to wall density, and the 


NACA TN 2110 


9 


boundary- layer thickness were chosen from experimentally measured 
values at an inlet pressure of 30 inches of mercury absolute and a 
total temperature of 560° E. A comparison of the assumed cosine 
density profile and the faired curve of the profile calculated from 
# the total-pressure measurements made at these conditions is shown 
in figure 9. The light path and error calculations for the assumed 
profile are shown in appendix B. From calculated light paths, the 
averaging, window, and displacement errors were each evaluated for 
a number of incident heights of light rays in order to obtain 
enough points to plot a distribution curve for each error. The 
sum of the three errors is the total error due to refraction. The 
four error distribution curves are plotted in figure 10. These 
results show that an uncorrected boundary- layer profile indicated 
by an interferometer for conditions such as the assumed ones is of 
little value in calculating aerodynamic data. By extending the 
calculations to other conditions of inlet pressure and temperature 
but still assuming the cosine density profile, the maximum total 
refraction error was found to be less than 1 percent in a tunnel 

span of 3.6 inches if the ratio /^Pi - P W A> is less than 0.19. 

Correction of refraction errors. - Analysis of boundary- layer- 
interference data is the reverse problem of the one previously con- 
sidered; that is, given an indicated density profile, find the 
density distribution through which the light passed. Unfortunately, 
the indicated profile provides insufficient information with which 
to make the necessary calculations. Any process of correction will 
therefore require certain assumptions and must be developed with the 
idea of producing a reasonably accurate result with a minimum of 
work. 


A possible method of correction that was considered makes use 
of the assumption that the indicated profile is a first approximation 
to the correct profile. Light paths calculated from a numerical 
integration of the refraction equation could be used to obtain values 
for the various errors that, when applied to the indicated profile, 
would give a corrected curve for the boundary- layer density distri- 
bution. Such a process would be of doubtful accuracy if the indi- 
cated profile was very different (5-percent maximum error in density) 
from the actual profile because of the sensitivity of the light path 
to values of the density gradient. For example, in figure 11, where 
the slope of the assumed and indicated profiles differ, the method 
is probably invalid. Also, the calculations would be lengthy unless 
programed for an automatic calculator. 
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For the case where time is not a factor in obtaining the inter- 4 

ference photographs, the following method, which was applied to the 
data presented, can be used: Two interference photographs were taken 

at each condition, one with the camera focused on the center of the 
tunnel and the other with the camera focused on the edge of the 
plate nearest the light source. The indicated profiles calculated 
from the two photographs differed only in that the apparent height 
above the surface at which each light ray appears to originate is 
shifted. The emergent angle 9 e of the light leaving the tunnel 
can be found at a given value of the density ratio by taking the 
difference in the apparent heights indicated for that density ratio 
by the two profiles and then dividing the difference by the dis- 
tance between the object planes for the two photographs. The angle 
9 e is related to the density difference between the incident and 
emergent points on the light path by the equation 



which was obtained by integrating equation (5), neglecting a vari- 
ation of the index of refraction in the z-direction (two-dimensional 
boundary layer). The two approximations made in deriving equation (7) 
are that tan 9 & = 9 Q and that n e 2 - n^ 2 = En^ng - n^). If the 

density is assumed to change linearly in the y-direction between the 
incident and emergent points on the light path, the following rela- 
tions that were found from the refraction equation for a constant 
density gradient (appendix C) can be applied: 



Substituting equations (7) and (9) in (8) yields 



NACA TN 2110 


11 


which is an expression for the averaging error. The window error 
can be found in terms of the fringe shift Ng from equation (4). 


A plot of the fringe shift N g against the angle 0 e of the light 

when it strikes the window for a l/2-inch thickness of crown glass 
(n g = 1.52) is shown in figure 12. The values of the two errors, 

when subtracted from the indicated density, give a number approxi- 
mately equal to the density ratio at the height y^. Equations (A8) 
(appendix A), (7), and (9) lead to an expression for y^ in terms 
of the apparent height and emergent angle, both of which are known 
from the photographic data: 


The fraction F, determined by the location of the object plane 
(fig. 8), is l/2 when y a is read from the curve corresponding to 
the center-focus interferogram. The resulting equation is obviously 
approximate because it only allows y.^ to be greater than or equal 


y^ can be less than y a . The approximation results because of the 
constant -density-gradient assumption. 

In order to determine the reliability of the assumptions involved 
in the derivation of the correction process, the process was applied 
to the indicated profiles calculated in appendix B. The agreement 
between the corrected points and the original assumed cosine density 
profile is a measure of the accuracy of the process . The corrected 
data points are shown in figure 13 with the assumed density profile 
and the indicated density profiles for the two locations of the 
object plane. 

The faired curve for the indicated profile with the camera 
focused on the edge of the plate (fig. 13) has been shifted to the 
left by 0.0023 inch from the curve shown in figure 11. The reason 
can be explained with the aid of figure 8. The intersection of the 
dashed line, which represents the backward extension of an emergent 
light ray, changes from the height y a at the center to a point 



' (ID 



( 12 ) 


to whereas the plot of 
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below the surface when, the object plane is shifted to the edge of 
the plate. Some of the light passing through the boundary layer 
therefore appears to come from a point below the surface of the 
plate. The distance below the surface depends on the emergent angle 
and the height of the ray. The calculations for the indicated den- 
sity profile in figure 11 showed that the lowest point from which 
the light would appear to originate would be 0.0023 inch below the 
shadow cast by the surface of the plate. In the evaluation of a 
photograph, the lowest point is considered as the surface, thereby 
introducing a displacement of the indicated profile equal to the 
shift of the surface. Because the y values used to calculate 
the emergent angle of the light are both referred to the same axis, 
the measured profile at the edge of the plate must be shifted to 
refer to the axis used in the evaluation of the indicated profile 
at the center of the tunnel. An interferogram showing the shift 
of the surface is shown in figure 14(a). The effect is markedly 
pronounced because a double exposure was made: one without flow, 

showing the actual location of the surface; and one with flow where 
the surface appears in the position to which it was shifted. With 
a laminar velocity profile, a double exposure is not required in 
order to detect the shift of the surface. Because of the zero den- 
sity gradient at the surface of the plate, the light that just 
grazes the surface will not be deviated in passing through the tun- 
nel. This light will therefore arrive at the same place on the 
negative both with and without flow. It shows up as a bright line 
or a discontinuity in the boundary- layer interference fringes. The 
discontinuity and the brightness just above the surface can be seen 
in figure 14(b). 

The corrected data points in figure 13 furnish insufficient 
information to plot a complete density profile. As a result of 
shifting the indicated profile in the aforementioned manner, values 
of the emergent angle of the light can only be found for the upper 
part of the boundary layer. Data with which to calculate the emer- 
gent angle in the region near the surface could be obtained by taking 
and evaluating a third photograph with the camera focused on the edge 
of the plate nearest the camera. Such a photograph could not be 
obtained in this investigation because of physical limitations of 
the equipment used. The behavior of the corrected density profile 
was therefore estimated in the region where data were unavailable. 

The results of an application of the correction process to an 
experimentally determined indicated density profile are shown in fig- 
ure 15, along with the indicated profiles for the camera focused on 
the center and edge of the plate. The inlet conditions under which 
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these data were obtained were a free- stream total pressure of 
29.33 inches of mercury absolute and a free- stream total tempera- 
ture of 560° B, which are very close to the conditions corre- 
sponding to the assumed profile. 

High-speed and normal- speed interference motion pictures were 
taken of the boundary layer in order to detect any fluctuations in 
the density profile, either instantaneously or over an extended 
period of time. The interference fringes of the boundary layer were 
steady in both cases for inlet pressures above 12.5 inches of mercury 
absolute. Below this pressure, the boundary layer separated from the 
surface and fluctuated up and down. This separation was caused by 
the starting shock, which had moved to the rear edge of the plate. 
Interference data have been evaluated only over the range where the 
boundary layer did not change with time. The process of taking two 
photographs at each condition over a finite interval of time in order 
to evaluate and to correct the interference data therefore did not 
cause any inaccuracies. A method has been suggested of adapting the 
correction process to the study of transient phenomena by using an 
infinite fringe adjustment of the interferometer and taking simul- 
taneous interferograms of two object planes. 

The calculations of the refraction errors and the correction 
process that is applied to the indicated density profile both make 
use of the assumption that light enters the boundary layer parallel 
to the surface. In actual practice, all the light cannot be parallel 
because the light source must have a finite extension. For the col- 
limating system used in this investigation, the maximum angle between 
any light ray and the central ray, which was adjusted parallel to the 
surface, was 0.00032 radian. In order to show that the light path 
through the boundary layer can be considered independent of the 
incident angle of the light, for small angles, the refraction equa- 
tion was integrated for an arbitrary density distribution. Differen- 
tiation of the resulting equation 



(13) 


with respect to 0^ shows that the percentage change of z with 
is of the same order of magnitude as 0^. The assumption of 
parallel light is therefore justified. 
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I nterferogram- evaluation errors, - Inaccuracies in measurements 
from an interferogram do not arise in the measuring instrument. They 
are dependent on the judgment of the operator in choosing the exact 
location of a particular point on the negative. The following effects 
add to this problem: 

1. In the region near the edge of the boundary layer, where the 
density gradient is relatively small, the diffuseness of the inter- 
ference fringes causes difficulty in locating points corresponding 
to fractional fringe shifts. The measurements in this region are 
therefore the least accurate of any of the distance measurements. 

2. The surface of the plate is not sharply defined in an inter- 
ferogram because of the angular deviation of the light from the 
extended source. A combination of reflected and defracted light 
from the plate interferes with the light from the beam that travels 
around the tunnel and the interference produces the pattern at the 
surface seen in figure 16. This region extends for a distance of 
0.0054 inch above the surface and causes an equivalent uncertainty 
in the location of the surface. In order to facilitate a more 
accurate estimate of the location of the surface, a photograph was 
taken of a known-size wire lying on the plate. By using the magni- 
fication of the picture and the wire size, the distance from the top 
of the image of the wire to the position at which the image of the 
surface should have been located was measured on the negative. 
Inasmuch as the appearance of the negative at this point was known, 
an estimate of the location of a similar point was made for each 
boundary-layer interferogram. The error is therefore much less 
than 0.0054 inch; the value depends on the judgmen J of the evaluator. 

3. Developing time for the interferograms should be so chosen 
that the bright and dark fringes are of equal width because meas- 
urement of fractional fringe shifts in evaluating the pictures is 
desirable. Overdeveloping shifts the boundary line between a bright 
and a dark fringe and widens the bright fringe, as shown in the 
interferogram of figure 17 . The location of the center of a fringe 
is unaffected by overdeveloping. 

4. Inasmuch as the interference measurements were made directly 
on the negative of the interferogram, no errors due to distortions 
in enlarging resulted. Distortions due to the optical system or to 
film shrinkage, however, are possible. By photographing a grid of 
uniform size and then measuring the distortion of the grid, the mag- 
nitudes of such errors were found to be negligible. The magnifica- 
tion of the photographs was also found from the photograph of the 
grid. 
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Tunnel -window boundary- layer errors. - The indicated density, 
calculated from equation (4) for a fringe shift N measured from 
an interf erogram, has been assumed to represent the sum of two 
terms according to the equation 


p ind 


av 



N. 


g 


(14) 


Because of the boundary layer on the tunnel windows, the assumption 
of a two-dimensional boundary layer on the plate used in arriving 
at equation (14) is not fulfilled. Equation (4) therefore yields 
a density 



J) 

kL 


(15) 


which differs from p^ n< ^. A 11 expression for this difference can be 

obtained by considering two rays of light, one traversing the free 
stream above the plate and the other passing through the boundary 
layer. If the same boundary layer on the windows is assumed for 
each light path, the average density traversed by each ray is 

for free stream, 


J S 


"\ 


P 1 (L - 26 g ) + 2 j p(z)dz 


p f = 

for boundary layer. 


0 


> (16) 


g 


P av (L - 26 g ) + 2 / p(z)dz 


P b = 


X 0 

+ — W K 

kL B J 


The fringe shift N used in equation (15) corresponds to the den- 
sity difference Pf ” Pb* Replacement of the last term of equa- 

tion (15) with the expression for p f - p b given by equation (16) 
yields 


- P 


av 


) f!s 

L 


p ind (^av + N g^ + ( p : 


(17) 
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Comparison of equations (14) and (17) shows that the density p'i n( fL 
used in the calculations is larger than p^ n(i hy the amount of the 

last term in equation (17), which introduces an equivalent error. 
Written in terms of density ratio, the error is 



Errors caused by end effects will therefore vary inversely with the 
span of the tunnel. The largest error will be for the ray traveling 
closest to the surface. For example, if the ratio of wall to stream 
density is 0.584 and the boundary layer on the tunnel window is 
l/lO inch thick, the largest error for a tunnel 3.6 inches wide will 
be 4.0 percent. The error in this calculation due to the assumption 
that the boundary layer on the windows is the same at all points 
will be a second- order effect. 

The tunnel width to be used for a particular boundary- layer 
investigation should be such as to make the sum of the refraction 
and the end-effect errors a minimum. Equation (8), which gives an 
order-of -magnitude estimate of the largest refraction error as being 
directly proportional to the square of the tunnel span, can be used 
to find the optimum width that balances the two errors when combined 
with equation (18). Calculations show that the tunnel width for this 
investigation should have been approximately 1.5 inches. Aerodynamic 
considerations may modify this calculation. 

Two-dimensional boundary layer. - Interferometer data will still 
be in error when corrected by the preceding methods if the boundary- 
layer flow is not two-dimensional. In order to determine how closely 
this condition of two- dimensionality was met in these experiments, 
total-pressure surveys were made of the boundary layer at three span- 
wise positions on the plate. Density profiles calculated from these 
surveys using equation (l), but assuming constant total temperature, are 
plotted in figure 18. The surveys were made at three Beynolds num- 
bers. The agreement between the three curves at each condition is an 
indication of the two-dimensionality of the boundary layer. Profiles 
for the lowest Beynolds number (fig. 18(a)) agree most closely. These 
profiles are probably an inaccurate representation of the actual den- 
sity distribution because of the assumption of constant total 
temperature. 


NACA TN 2110 


17 


Velocity-profile calculations. - Velocity profiles were cal- 
culated from the interference data by assuming a theoretical total- 
temperature distribution, as was done for the total-pressure data. 

The calculation differs slightly in one respect. An accurate value 
of density at the wall could not be obtained from an interferogram. 

The wall density was therefore taken to be the value corresponding 
to the measured surface temperature. According to the energy equa- 
tion, however, surface temperature is determined by the local stream 
Mach number M-^ and the temperature recovery factor by the 

equation 

£ - 1 - °- 2 ( 19 ) 

X 1 1 + 0.2 M -yf 

A value of of 0.845 for a Prandtl number Pr of 0.72 is found 

in reference 6. Inasmuch as the measured value of surface tempera- 
ture indicated a Prandtl number higher than 0.72, the wall tempera- 
ture calculated from the measured M]_ with equation (19), using 

t)-j- = 0.845, was lower than the measured surface temperature. From 

equations (l) and (3), the condition for zero velocity at the wall 
is 

~~ (1 + 0.2 M - l 2 ) = 1 ( 20 ) 

T 1 Pi 

ft 

In order to use the curves given in reference 6 and at the same time 
satisfy equation (20) for the measured values of P v /p-± and M^, 

it was therefore necessary to adjust the value of M]_ used in equa- 
tion (19). The new value of M^ was used to calculate the total- 

temperature distribution from the theory. This value was approximately 
0.05 larger than the measured M^. The measured values of p/p^ 
and M^ and the calculated total- temperature distribution furnish 

the necessary quantities to calculate Mach number distribution. The 
velocity ratio was calculated from the density and the Mach number 
curves using equations (2) and (3). 


RESULTS AND DISCUSSION 
Velocity and Density Profiles 

The data needed to calculate velocity and density profiles were 
obtained by interferometer and pressure-probe measurements in the 
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local Reynolds number range Re x of 0.318 X 10® to 1.08 X 10®. 

Mach number changed slightly over this range of operating condi- 
tions; the nominal value was 2.02. The interference photographs 
taken with the camera focused on the center of the tunnel are shown 
in figure 19. Indicated density profiles measured from these photo- 


graphs at a position 2^ inches from the leading edge of the plate 

u 

and the corresponding corrected profiles are plotted in figures 20(a) 
and 20(b), respectively. The theoretically more interesting plot of 
the ratio of density in the boundary layer to density at the edge of 
the boundary layer against the dimensionless distance variable 



is shown in figure 21(a) for the interference data and in 


figure 21(b) for the total-pressure data. These figures include a 
curve calculated from the theory presented in reference 6 for the 
case of an unheated, uncooled plate in a supersonic stream of Mach 
number 2. Similar experimental and theoretical curves for the 
velocity profiles are shown in figures 22(a) and 22(b). 


The density profiles calculated from the total-pressure data 
are nearly the same in shape as the theoretically calculated curve. 
The shifting of the experimental profiles to the right of the theo- 
retical curve might be due to a random error in the experiments 
because no apparent relation exists between the shift of the profiles 
and the variables pressure, temperature, or Reynolds number. 


Relatively good agreement exists between theory and the profiles 
obtained from pressure measurements. The interferometer density 
profiles are below the pressure measurements and theory near the 
edge of the boundary layer. The lower density indicated for this 
region may be connected with the fact that the interferometer meas- 
ures the thermal boundary layer, whereas the total-pressure probe 
measures the velocity layer. The increasing inaccuracy of the indi- 
cated profiles and of the correction process at the higher inlet 
pressures, however, is more probably the explanation of the dis- 
crepancy. A comparison of figures 22(a) and 22(b) shows that the 
profile for the lowest inlet pressure is in close agreement with 
the total-pressure profile. The agreement of the other curves 
decreases as the inlet pressure increases. Application of the cor- 
rection process to the indicated profile calculated for the assumed 
density variation produced a corrected curve having a maximum error 
of 1.2 percent. All the factors that influence the indicated dens- 
ity profile in actual practice have probably not been accounted for, 
inasmuch as the errors in the profiles of figure 21(a) are apparently 
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greater than this small percentage. For example, the effect of 
boundary layer on the tunnel windows is disregarded in both the 
evaluation and the correction of the indicated profile. 

The boundary- layer thickness taken from the interference data 
is considerably greater than that indicated by the pressure data. 

The reason is apparent when the boundary layer in the corners between 
the plate and the windows is considered, as in reference 7. The bound- 
ary layer on the plate and on the window must form a continuous var- 
iation of the density around the corner rather than the abrupt right- 
angle change that would exist if the two boundary layers did not 
interact with each other. This interaction has the effect of thick- 
ening the boundary layer on the plate near the corners. A ray of 
light that enters the tunnel at a point just at the edge of the 
boundary layer will therefore travel a shorter optical path than a 
ray passing through at a height high enough to be unaffected by the 
corner boundary layer because of the reduced density in the corner 
regions that extends above the edge of the boundary layer . Inter- 
ference measurements are therefore useful only for an approximate 
estimate of the boundary- layer thickness. 

Boundary- layer thickness taken from the density profiles of the 
pressure measurements is plotted in figure 23 with a curve calculated 
from the theory of reference 6. The theoretical curve gives values 
approximately 30 percent lower than the measured values. The exper- 
imental points, however, are insufficient to establish a trend. 


Skin Friction 

The velocity and density profiles were used to calculate the 

section skin-friction coefficient C f to a position 2- inches from 

2 

^he leading edge of the plate by a graphical integration of the 
Karman momentum equation, which was used in the following form: 


C f 


1 2 - 
2 P 1 U 1 ^ 


r\b 


2 

x 


_P_ JL 
P 1 "i 



( 21 ) 


Values of are plotted in figure 24 as a function of local stream 

Reynolds number. These points are compared with curves from the 
incompressible and compressible laminar -boundary- layer theory. The 
results indicate that the laminar skin-friction coefficient increases 
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with Mach number, rather than decreasing as predicted by theory. 

The experimental conditions are inconsistent with the assumptions 
of the theory in that a disturbance produces rotational flow at the 
leading edge. This rotational flow may be equivalent to starting 
with an initial boundary layer at the leading edge. If such is 
true, the value of x used to calculate Re x and should be 


larger than the actual distance (2^ in.) by an amount that would be 

u 


= 3— inches. 
2 


required to develop the initial boundary layer without the leading- 
edge shock. As an illustration, the point Re x = 0.61 X 10®, 

Cf = 0.00205 (fig. 24) has been recalculated for x 

The new values are Re x = 0.72 X 10^ and Cf = 0.00147, which is 
a point slightly below the compressible-theory curve. The meas- 
ured values of Cf agree with the theory in that they decrease with 

The agreement between the values of the skin- 


an increase in 


A/Ei 


x # 


friction coefficient calculated from the total-pressure and inter- 
ferometer measurements , if used as a measure of accuracy of the 
interference method, shows that the correction process overcorrects 
when light refraction is small (low Reynolds number, fig. 24) and 
undercorrects when the refraction is large. The interference meas- 
urements give a value that is 9 percent below the total -pressure 
measurements at the lowest Reynolds number and 5 percent above at 
the highest Reynolds number. 


Recovery Factor 

The recovery factor r\^ f which is the fraction of the kinetic 
energy of the undisturbed flow that raises the effective temperature 
of the plate to a value above the ambient temperature, was calculated 
from the relation 


t-v - ti 



( 22 ) 


which was obtained from the definition of total temperature and the 
equation defining t)^. The temperatures t v and Tq were measured 

and t^ was computed from M]_ and Tq. Although the recovery factor 

is theoretically dependent only on the Prandtl number of the air, the 
results of these measurements showed a tendency for to increase 

with Reynolds number. The experimentally determined values are 
plotted in figure 25 as a function of the Reynolds number in the test 
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section. An integration of the energy equation for incompressible 
laminar flat-plate flow, performed by Pohlhausen (reference 8), 
predicts a recovery factor of 0.849, which is the square root of 
the Prandtl number of 0.72. (Reference 6 gives 0.845 for Prandtl 
number of 0 . 72 .) The measured values are slightly larger than 0.87 
but not as large as 0.90, which has been found for turbulent flow. 

A recovery factor of 0.87 corresponds to a Prandtl number of 0.76. 


Errors due to light refraction must be considered in any boundary- 
layer investigation using light- interference measurements. A laminar 
boundary layer on an unheated, uncooled plate has been shown to cause 
a refraction error of less than 1 percent for a tunnel 3.6 inches 


the edge of the boundary layer, P w is the density at the surface 

of the plate, and 5 is the boundary -layer thickness. Theoretically, 
density profiles calculated from interference data can be corrected 
for refraction effects so as to have a maximum error in the boundary - 
layer density of a few percent. A comparison between values of the 
skin-friction coefficient calculated from pressure and from corrected 
interference measurements shows that the error is generally larger 
than the theoretically indicated maximum. Because refraction errors 
are directly proportional to the square of the tunnel span, whereas 
errors due to end effects are inversely proportional to the span, 
an optimum tunnel width exists for which the total error caused by 
both these errors will be a minimum. The tunnel used for this 
investigation was larger than optimum width. 

Density and velocity profiles calculated from total -pressure 
data were similar in shape to profiles predicted by laminar -boundary - 
layer theory. Profiles from interference data were in poor agreement 
with the theory near the edge of the boundary layer. The two methods 
of measurement had the best mutual agreement at the lowest Reynolds 
numbers used in the investigation. Total -pressure -tube measurements 
gave the most reliable value of boundary -layer thickness. Skin- 
friction coefficients, calculated with the momentum equation, varied 
inversely with the square root of Reynolds number. The values of 
the section skin -friction coefficient for a given Reynolds number 
were higher than for incompressible flow, rather than lower as pre- 
dicted by theory. Consideration of the leading-edge shock wave 
offers a possible explanation of this discrepancy . The recovery 
factor calculated from measured temperature and Mach number showed 
a tendency to increase with Reynolds number. The values were in the 


CONCLUSIONS 


wide if the ratio „ is less than 0.19, where is the density at 
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region of 0.87 to 0.88, which was slightly larger than the value 
calculated from an integration of the energy equation for a laminar 
incompressible boundary layer. 


Lewis Flight Propulsion Laboratory, 

National Advisory Committee for Aeronautics, 
Cleveland, Ohio, November 23, 1949. 
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APPENDIX A 
SYMBOLS 

The following symbols are used in this report: 

Cf section skin- friction coefficient 

Df skin-friction drag force on plate section 

F fractional portion of tunnel span between object plane and 
tunnel window (fig. 8) 

k constant (k = 0.1166 cu ft/slug for air when = 5461 A) 

L span of tunnel 

M Mach number 

N fringe shift through boundary layer measured from interferogram 

Ng fringe shift caused by light striking tunnel window at angle 

n refractive index of air 

ng refractive index of tunnel window 

P total pressure 

Pr Prandtl number 

p static pressure 

Re x local stream Reynolds number based on distance x 

S length along light path through boundary layer 

T total temperature 

t ambient-air temperature 

u air velocity 

chordwise distance measured from leading edge parallel to 
stream 


x 
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y 

distance measured perpendicular to plate surface 

y' 

dy/dz 

y" 

d 2 y/dz 2 

z 

transverse coordinate between tunnel windows 

a 

nondimens ional coordinate perpendicular to plate surface 

3 

nondimens ional coordinate parallel to plate surface and 
transverse to stream direction 

6 

boundary- layer thickness 

6 g 

boundary- layer thickness on tunnel windows 


temperature recovery factor 

9 

angle between tangent to light path and plate surface 

d v 

refraction angle of light impinging on tunnel window 


vacuum wavelength of light 


stream kinematic viscosity 

P 

air density 

(p/p^ ) av 

density ratio averaged along light path between tunnel 
windows for two- dimens ional-flow field 

(p/oi)ind 

density ratio calculated from fringe shift N with 
equation (4) 

T 

thickness of tunnel window 


Subscripts : 


0 

free- stream stagnation conditions 

1 

conditions at edge of boundary layer 


a 


conditions at point where light ray appears to originate 
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av 

b 

e 

f 

i 

ind 

w 


average 

average density along light path through boundary layer 
including side-wall boundary layer 

condition at point where light path intersects tunnel 
window nearest camera 

average density along light path through free stream 
including side-wall boundary layer 

conditions at point where light ray is incident on 
boundary layer 

indicated 

conditions at surface of plate 
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APPENDIX B 


LIGHT PATH AND ERROB CALCULATIONS FOB 


ASSUMED PROFILE 


Equations (5) can be put in the form 


d£ _ 1 dn 
dz ~ n dy 


(Bl) 


if n is independent of z; 9 is the angle that a tangent to a 
point on the light path makes with the z-axis. By substituting 
density for refractive index, according to equation (6), and replac- 
ing z and y with the nondimens ional quantities a = y/6 and 
S = z/L, where 6 is the boundary- layer thickness and L the 
tunnel span, equation (Bl) becomes 


a 2 » 


i a (£) 


dp‘ 


da 


(B2) 


Because 6 is small, it has been replaced by dy/dz and l/n has 
been assumed 1. The assumed density profile is represented by 


_P _ 1 

Pi -2 




cos ita 


(B3) 


where p w is the density at the wall and 


dl t> . 


\ p l/ it f 1 p w . 

■fcT ' 2 - — I «m «a 


(B4) 


Equation (B4) substituted into equation (B2) gives the differential 
equation for the light path. 


d 2 a B 2 kjt(pi~ p w ) 


sin ita = K sin na 


(B5) 


dp 
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When integrated once, equation (B5) yields 


da , l~2K , ,1/2 

^ = /W— («c 1 - cos ita) (B6) 

If the light enters the tunnel parallel to the surface of the plate, 
<Iz = £ ^ = ® when a = a^, the initial height of the light on 
entering the tunnel. By replacing the constant of integration c-^ 
in equation (B6) and integrating again, 

r\a 

p = vs 


a- 


4 


da 


cos jra^ - cos Tta 


(B7) 


where 


a > 


This equation is an elliptic integral of the first kind, which has 
been evaluated for a number of values of and is plotted in 

figure 26 in terms of a and a/jtKP. From this information, a e 
and 0 e , the height above the surface and the angle of the light 

when it leaves the tunnel, respectively, were obtained as functions 
of for a free-stream total pressure of 30 inches mercury abso- 

lute, and a free-stream total temperature of 560° R; that is, 

= 5.15 X 10“ 5 slugs per cubic foot, P w /pp = 0.584, 

6 = 0.028 inch, and L = 3.60 inches. Curves showing a g and 0 e 
are presented in figure 27. The apparent height at which each 
light ray appears to originate was calculated from and a e 

for two locations of the object plane, one on the center line of 
the plate, the other at the edge of the plate nearest the light 
source, from the equation 



n 


g 


0e 


y a = 6 a, e 


FL + T 


(B8) 
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where F is determined by the location of the object plane. The 
indicated density was found by numerically summing the density along 
each light path to obtain the average density traversed by the ray 


and then adding the term 


*0 

PpkL 


N 


g 


due to the fringe shift when the 


light passes through the tunnel window at an angle. The indicated 
profiles focused on the two planes and the original assumed profile 
are shown in figure 11. 


The manner in which the error between the indicated and assumed 
profiles varies through the boundary layer is plotted in figure 10. 
The total error is composed of the three errors previously discussed 
in Befraction errors. Each of these errors was evaluated from the 
light-path calculations in order to show the individual contribu- 
tions to the total error. 


The averaging error was found by subtracting the numerically 
calculated average density along a light path from the density at 
the incident height of the light. Thus the averaging error is 



(B9) 


higher than 


The window error caused by the light striking the tunnel win- 
dow at an angle makes the indicated density ( ■£- \ 

\ p l/ind 

the average density. The measured fringe shift N, used to cal- 
culate the indicated density, will be less by an amount Ng than 


the fringe shift required to indicate a density 




g 

Therefore, 



= 1 - 


av 




or 



ind 



\> 

PnkL 


av 


(ii) 
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i 


which can he used to calculate the window error by using the appro- 
priate value of Ng corresponding to the emergent angle of the 

light. 

The displacement error is given by the difference between the 
densities at the incident height and the height y a from which 

the light appears to originate. The sign of the error 



depends on whether- y a is greater or less than y^. 

Summing the expressions for the three errors yields the identity 



(BIO) 


30 


NACA TN 2110 


APPENDIX C 


DERIVATION OF AVERAGE- DENSITY AND EMERGENT-HEIGHT 


EQUATIONS FOR CONSTANT DENSITY GRADIENT 

Equation (5) can be put in the form 

d£ = 1 dn 
dz n dy 


(Bl) 


Rewriting this expression in terms of density by means of equa- 
tion (6), including the assumption that l/n =1, leads to a dif- 
ferential equation that can be integrated for a constant density 
gradient, giving 

6 = k z (Cl) 

dy 

where the boundary conditions at the point of entrance of the light 

ray are 0=0° and z = 0. Because 0 will be small, 

tan 0=0= dy/dz. Equation (Cl), when integrated again, becomes 


V - Ji = 



Z 2 


(C2) 


The light ray enters the tunnel at the height y ^ at z = 0 and 
leaves at y e at z = L. Under the assumption of a linear density 
change from y e to yi, the density gradient can be expressed as 


dp = 
dy 



(C3 ) 


Substituting equation (C3) in equation (C2) and evaluating at 
z = L gives 


(y e - Yi) 


2 




( 9 ) 
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The average density along a light path is 


p av - | 


p(S) dS 

0 


(C4) 


Inasmuch as dS = dz 1 + y ' ^ = dz 'sj 1 + 9^ , and 9 ^ is negli- 
gible compared to -unity, dS can be replaced vith dz. The assump- 
tion of constant density gradient is expressed in the following 
equation: 


p(y) = p w + ^ y (C5) 

When rewritten as a function of z using equation (C2 ) , equa- 
tion (C5) becomes 


p(z) = p w + § y i 



(C6) 


(Note that p w + ^ y^ is equivalent to p^) Substituting this 

expression for density variation along the light path in equa- 
tion (C4) and integrating yields 



(C7) 


which can be rearranged into equation (8) using equation (C3). 
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Figure 1. - Interferometer and wind-tunnel installation. 
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Figure 2. - Wind tunnel and flat plate with one side plate removed. 
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Figure 3. - Typical opening of boundary-layer total -pressure probe. 
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Figure 4. - Leading-edge shock wave. 
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Distance from leading edge, x, in. 

Figure 5. - Pressure distribution on plate at various inlet 
pressures and free-stream total temperature of 560° R. 


Total-temperature ratio, T/Tq 
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Figure 6. - Boundary -layer total-temperature distribution calculation from theory of 
reference 6. Mach number at edge of boundary layer, 2.09. 
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Fringe shift 



Figure 7. - Illustration of method for analysis of interference data. 
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Figure 8. - Schematic diagram of light path through boundary layer. 


Ol 


NACA TN 21 10 


Density ratio, p/p 



Figure 9* - Comparison of measured and assumed density profile used for refraction 

calculations. Free-stream total pressure, 30.00 inches mercury absolute; free-stream 
total temperature, 560° R; Mach number at boundary -layer edge, 2.01; local stream 
Reynolds number, 0.759 x 10$. 
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Figure 10. - Variation of refraction errors through assumed boundary -layer profile 
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Figure 11. - Comparison of indicated and assumed density profiles 
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Fringe shift caused hy tunnel window 
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0 12 24 36 48 60 


Incident angle of light at window, 0 e , min 

Figure 12. - Fringe shift caused by light striking l/2-inch 

thickness of crown glass (n = 1.52) at an angle. 
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Density ratio, p/p 



Figure 13. - Comparison of corrected interference data with assumed density profile. 
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(a) Double exposure. 



(b) Single exposure. 

Figure 14. - Interferograms illustrating apparent dovnward shift of surface when camera 

is focused on edge of plate. 




Density ratio, p/p 



Figure 15. - Comparison of indicated and corrected density profiles. Free-stream total 
pressure, 29.33 inches mercury absolute; free-stream total temperature, 5600 r- Mach 
number at boundary -layer edge, 2.04; local stream Reynolds number, 0.730 x 10^. S 
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Figure 16. - No -flow Interference photograph showing poor definition of plate surface. 



Figure 17 . - Interference photograph illustrating unequal thickness of dark and light 

fringes caused by overdeveloping . 
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Density ratio, p/p 



(a) Reynolds number, 0.388 x 10 6 . 


Figure 18. - Comparison of density profiles at three spanwise positions. 
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Density ratio, p/p 
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(b ) Reynolds number, 0.783 x 10 6 . 
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Figure 18. - Continued. Comparison of density profiles at three spanwise positions. 
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(c) Reynolds number, 1.09 x 10 6 . 




Comparison of density profiles at three 


spanwise positions. 
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Figure 18. - Concluded 
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(a) Free -stream total pressure, 15.93 inches mercury; stream Mach number 2.06; Reynolds 

number, 0.384 x 10°. 



(b) Free -stream total pressure, 19.52 inches mercury; stream Mach number, 2.05; Reynolds 

number, 0.483 x 10 6 . 



(c) Free -stream total pressure, 24.71 inches mercury; stream Mach number, 2.05; Reynolds 

number, 0.611 x 10°. 

Figure 19. - Interference photographs for six Reynolds numbers taken with camera focused 
on tunnel center line. Free -stream total temperature, 560° R. 
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(d) Free -stream total pressure, 29.33 inches mercury; stream Mach number, 2.04; Reynolds 

number, 0.730 x 10 6 . 



(e) Free -stream total pressure, 34.41 inches mercury; stream Mach number, 2.03, Reynolds 

number, 0.858 x 10 6 . 



(f) Free-stream total pressure, 39.55 inches mercury; stream Mach number, 2.02; Reynolds 

number, 0.989 x 10°. 

Figure 19. - Concluded. Interference photographs for six Reynolds numbers taken vith 
camera focused on tunnel center line. Free-stream total temperature, 560° R. 
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(a) Indicated profiles 2-g Inches from leading edge. Taken from interferograms of figure 19. 


2: 

> 

o 

> 


ro 

o 


Figure 20. - Density profiles. 
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Figure 21. - Density profiles from Interference and total-pressure measurements. 
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. Density profiles from interference and total -pres sure measurements. 


Figure 21. - Concluded 
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Figure 22. - Velocity profiles from interference and total-pressure measurements. 
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Figure 22. - Concluded. Velocity profiles from interference and total-pressure 

measurements. 
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Figure 23. - Thermal boundary -layer thickness obtained from total-pressure measurements 
as function of local stream Reynolds number; x = 2^ inches. 
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Figure 24. - Section skin-friction coefficient as function of local stream Reynolds 
number calculated from interference and total-pressure measurements; x = 2 ^ inches. 
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. - Temperature recovery factor as function of Reynolds 
number in test section. 
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figure 26. - a as function of 'fix. p 

VSic p = £ /’ a — — 

A/2 Vcos rra^ - cos Tra 


evaluated from equation 
da 


(B7). 
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Figure 27. - Emergent height and angle of light as function of entrance 

height. 
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